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ABSTRACT 

We have developed and characterized an imaging instrument to measure the spatial properties of 
the diffuse near-infrared extragalactic background light in a search for fluctuations from z > 6 galaxies 
during the epoch of reionization. The instrument is part of the Cosmic Infrared Background Experi- 
ment (CIBER), designed to observe the extragalactic background light above the Earth's atmosphere 
during a suborbital sounding rocket flight. The imaging instrument incorporates a 2° x 2° field of view, 
to measure fluctuations over the predicted peak of the spatial power spectrum at 10 arcminutes, and 
7" x 7" pixels, to remove lower redshift galaxies to a depth sufficient to reduce the low-rcdshift galaxy 
clustering foreground below instrumental sensitivity. The imaging instrument employs two cameras 
with AA/A ~ 0.5 bandpasscs centered at 1.1 ^m and 1.6 /im to spectrally discriminate reionization 
extragalactic background fluctuations from local foreground fluctuations. CIBER operates at wave- 
lengths where the electromagnetic spectrum of the reionization extragalactic background is thought 
to peak, and complements fluctuations measurements by AKARI and Spitzer at longer wavelengths. 
We have characterized the instrument in the laboratory, including measurements of the sensitivity, 
flat-field response, stray light performance, and noise properties. Several modifications were made to 
the instrument following a first flight in 2009 February. The instrument performed to specifications 
in subsequent flights in 2010 July and 2012 March, and the scientific data are now being analyzed. 

Subject headings: (Cosmology:) dark ages, reionization, first stars - (Cosmology:) diffuse radiation 
- Infrared: diffuse background - Instrumentation: photometers - space vehicles: 
instruments 



1. INTRODUCTION 

The extragalactic background light (EBL) is a measure 
of the integrated radiation produced by stellar nucleosyn- 
thesis and gravitational accretion over cosmic history. 
The EBL must contain the radiation produced during 
the epoch of reionization (the reionization-EBL, or sim- 
ply the REBL). The REBL comes from the UV and op- 
tical photons emitted by the first ionizing stars and stel- 
lar remnants, radiation that is now rcdshifted into the 
near-infrared (NIR). The REBL is expected to peak at 
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1-2 um due to the redshiftcd Lyman-a and Lyman-break 
features. Furthermore while the brightness of the REBL 
must be sufficient to initiate and sustai n ionization, the 
indiv idual sources may be quite faint (|Salvaterra et all 

i2oTa . 

We have developed a specialized imaging instrument 
to measure REBL spatial fluctuations, consisting of two 
wide-field cameras that arc part o f the Cosmic Infr ared 
Background Experiment fCIBER; iBock et al.ll2006h . de- 
veloped to measure the absolute intensity, spectrum, and 
spatial properties of the EBL. CIBER's imaging cameras 
are combined with a low-resolution spectrometer (LRS; 
iTsumura et alJl2~012t ) designed to measure the absolute 
sky brightness at wavelengths 0.75 < A < 2.1 /xm, and 
a narrow-band spectrometer (NBS; iKorngut et al.ll2~012D 
designed to measure the absolute ZL intensity using the 
854.2 nm Call Fraunhofcr line. A full description of the 
CIBER payload, including the overall mechanical and 
thermal design, and detailed descriptions of the focal 
plane housings, calibration lamps, shutters, electronic 
systems, telemetry and data handling, laboratory cali- 
bration equipment, flight events, and fligh t thermal per- 
formance, is given in IZemcov et all (|2012h . The observa- 
tion sequen ce and science targets from the first flight are 
available in ITsumura et all (|2010f ) . 

In this paper, we describe the scientific background of 
EBL fluctuation measurements in sections 11.11 and 11.21 
the instrument design in section^ laboratory instrument 
characterization in section [3l modifications following the 
first flight in section [4j and performance in the second 
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flight in section [SJ Sensitivity calculations are given in a 
short appendix. 

1.1. Science Background 

Searching for the REBL appears to be more tractable 
in a multi-color fluctuations measurement than by ab- 
solute photometry. Absolute photometry, measuring 
the sky brightness with a photometer and remov- 
ing local foregrounds, has proven to be problematic 
in the NIR, where the main difficulty is subtract- 
ing the Zodiacal light (ZL) foreground, which is a 
combination of scattered sunlight and thermal emis- 
sion from interplanetary dust grains in our solar sys- 
tem. However, absolute photomet ry studies give con- 
sistent resul t s in t he far-infrared faauser etall H998l 
Fixsen et all I1998L I Juvela et all 120091 iMatsuura et al l 
201 1L IP enin et aLl 12011ft . These far-infrared mea- 



surements are close to the EBL derived from galaxy 
counts though statistical and len sing techniques that 
probe below the confusion limit ( Mars den et al.l 120091 
Zemcov et all 120101 iBethermin et all 120101 iBerta et all 



2010ft" '. However in the NIR, at wavelengths appropri- 
ate for a REBL search, absolu te EBL measurements 
are not internally c o nsistent (|Cambresv et al.1 I2001L 
Dwek fc Arendtl 119981 IMatsumoto et al.l 120051 IWrightl 
20011 iLevenson fc Wright! 12008ft . A significant compo- 



nent of this disagreement i s related to the choice of 
mode l used to subtract ZL (jKelsall et al.l 119981 IWrightl 
12001ft . Furthermore, some absolute EBL measure - 
mcnts (jCambresv et all I2001L IMatsumoto et al.l [2005) 
are significantly higher than t he integrated galaxy light 
derived from source counts (iMadau fc Pozzettil I2000I. 
Totani et all 120011 ILevenson et al.l 120071 iKeenan et al.l 



2010ft~ 



The current disagreement between absolute measure- 
ments and galaxy counts are difficult to recon cile with 
theoretical calculations ([Madau fc Silkl 12005ft or TeV 



absorption measurements from blazars ( Gilmorel I2001L 
lAharonian et al.l 120061 ISchroedterl 12005ft . However TeV 
constraints on the NIR EBL req uire an assumption about 
the intrinsic blazar spectrum (jDwek et al.1 12005ft . Fur- 
thermore cosmic rays produced at the blazar arc not 
attenuated by the EBL and can produce secondary 
gamma rays that may explain the current TeV data 
without placing a serio us constraint on the NIR EBL 
(|Essev fc Kusenkd[2OT0ft . 

Instead of measuring the absolute sky brightness, it is 
possible to detect or constrain the REB L by studying 
the spatial properties of the background dCoorav et al.l 
120041 iKashlinskv et al.ll200~4ft . A spatial power spectrum 
of the EBL contains a REBL clustering component, evi- 
dent at an angular scale of approximately 10 arcminutcs 
as shown in Figure Q] that is related to the underlying 
power spectrum of dark matter. Numerical simulations 
of first galaxy formation i ndicate the effects of no n-linear 
clustering are significant ([Fernandez et al.ll2010l ). There 
are also REBL fluctuations from the Poisson (unclus- 
tered shot noise) component, but the amplitude of this 
term is more difficult to predict as it is related to the 
number counts of the first galaxies, that is, the bright- 
ness distribution and surface density of sources. In addi- 
tion, REBL fluctuations are thought to have a character- 
istic electromagnetic spectrum, peaking at the redshift- 
intcgratcd Lyman-a emission feature. If reionization oc- 



curs at z ~ 10, this emission peak is redshifted into the 
NIR, with a spectral shape that depends on the luminos- 
ity and duration of the epoch of reionization. 

Early measurements with th e Diffuse Infrared Back- 
groun d Experiment (DIRBE; Kashlins kv fc Odenwaldl 
2000ft and the Infrar ed Telescope in Space (IRTS; 



Matsumoto et all 12005ft used fluctuations as a tracer of 



the total EB L. A first detection of RE BL fluctuations was 
reported by IKashlinskv et al.l ( 2005ft using the Spitzer 
Infrared Array Camera (IRAC: [Fazio et al 112004ft in the 
3.6 and 4.5 /im bands in 5 x 5 arcminute regions, cor- 
responding to the IRAC field of view. The authors ob- 
serve a departure from Poisson noise on 1—5 arcminute 
scales which they attribute to first-light galaxies, after 
ruling out Zodiacal, Galactic, and galaxy clustering fore- 
grounds. The observed brightness of the fluctuations is 
approximately constant at 3.6 and 4.5 fim. This anal- 
ysis was later extended to 10 x lOar c min fields, giving 
similar results (Kashlin skv et al1l2007l ). I Thompson et al.l 
(|2007aft studied a 144 x 144 arcsec field with the Hubble 
Space Telescope (HST) at 1.1 and 1.6 /im, finding no ev- 
idence for z > 8 galaxie s contributing to the H ST or 
the Spitzer fluctuations ([Thompson et all l2007bft . Fi- 
nally, [Matsumotoler^l] ([2011ft report first-light galaxy 
fluctuations with AKARI at 2.4, 3.2 and 4.1 fim in a 10 
arcminute field. Their reported spectrum shows a strong 
increase from 4.1 to 2.4/xm, consistent with a Rayleigh- 
Jcans spectrum. 

In Figure Q] we show two predictions related to the 
angular power spectrum of REBL aniso tropics. The 
lower prediction (solid red line) is from iCoorav et al.l 
([2012ft . derived from the observed luminosity functions 
of Lyman dropout g alaxies at redshifts of 6, 7 and 8 
([Bouwens et al.1 [2008ft at the bright end. The reioniza- 
tion history involves an optical depth to electron scat- 
tering of 0.09, consistent with th e WMAP 7- year mea - 
surement of r = 0.088 ± 0.014 dKomatsu et'alJ 12011ft . 
The absolute REBL background is 0.3nW m~ 2 sr~ 1 at 
3.6 /xm for this model . ICoorav et al.l (120121 ) improved on 
previous predictions ([Coorav et al.l 12004ft by accounting 
for non-linear clustering at small angular scales with a 
halo model for reionization galaxies at z > 6. Note that 
the REBL fluctuation power is similar at 1.6 and 1.1 /im 
given the redshift of reionization is around at z ~ 10. 

The upper prediction (dashed red line) is normalized to 
the anisot ropy amplitude level rep orted by Spitzer-IRAC 
at 3.6 /im (IKashlinskv et al.ll2005ft . This power spectrum 
requires an absolute REBL background between 2 to 3 
nWm~ 2 sr~ 1 at 3.6 /xm. We scale the power spectra to 
shorter wavelengths based on a Rayleigh- Jeans spectrum, 
consistent wi th the combined measu rements of Spitzer 
and AKARI (|Matsumoto et al.|[20Tlft . 

Fluctuation measurements are only feasible if the 
contributions from foregrounds can be removed. For- 
tunately, it appears easier to remove foregrounds in 
fluctuation measurements than in absolute photometry 
measurements. The largest foreground, ZL, is known 
to be spa tially uniform on spatial scales smaller than 
a degree ([Abraham et al.l 119971 IKashlinskv et al.l 120051 
IPvo fc et al.l 120 lift . Furthermore, any spatial variations 
in ZL can be monitored and removed by observing a field 
over a period of time, as the view through the interplane- 
tary dust cloud changes annually. Galaxies and stars give 
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Fig. 1. — Power spectra of REBL and foreground fluctuations at 1.1 (im (left) and 1.6 fim (right). In both cases the clustering power 
spectra of local (z < 3) galaxies, for sources brighter than two different magnitude cutoffs, arc shown as the bl ue solid and dashed lines. 
These galaxy clustering power spectra are based on measured fluctuations as a function of cutoff magnitude from Sullivan et al. (2007) and 
are consistent with the predictions by Hcleason ct al. (2012) based on a large compilation of galaxy luminosity functio ns b etween z = 
and 4. The two red lines correspond to two expectations on the REBL anisotropy power spectrum as described in scction ri.il Upper limits 
to the ZL fluctuation power, shown in black, are scaled from experimental upper limits at longer wavelengths by the ZL spectrum. The 
predicted CIBER sensitivities in both bands are shown in orange. These are calculated with the instrument parameters listed in Table 
[2] assu ming that the detector noise given in Table [4] is uncorrelated and Gaussian over the array and using the AC^ formalism in IKnoxl 



spatial fluctuations from Poisson variations and cluster- 
ing. These can be eliminated by masking sources from 
the image, either through detection or by using an ex- 
ternal catalog of known sources. Galaxy clustering, ar- 
guably the most serious of these potential contaminants, 
requires a sufficiently deep source cutoff to reduce the 
clustering spectrum below the level of REBL fluctuations 
by masking sources. 

1.2. Theoretical Design Drivers 

These early fluctuation results call for a next genera- 
tion of improved measurements at shorter wavelengths, 
spanning the expected peak of the REBL electromag- 
netic spectrum, with wide angular coverage, to defini- 
tively measure the expected peak in the REBL spatial 
power spectrum. In order to make a definitive REBL 
fluctuations measurement, we require: (1) a wide field of 
view to allow measurements of the characteristic REBL 
spatial power spectrum, (2) observations in multiple NIR 
bands in order to characterize the REBL electromagnetic 
spectrum and distinguish it from potential foregrounds, 
and (3) arcsecond angular resolution to remove galaxies 
to a sufficient depth to minimize the galaxy clustering 
foreground signal. 

High-fidelity spatial imaging on degree scales is prob- 
lematic in the NIR due to airglow emission from the 
Earth's atmosphere, which is some 200—1500 times 
brighter than the astroph ysical sky in the NIR J, H 
and K bands (IAllenlll976f). Airglow e mission has time- 
variable structure (jRamsav et al.lll992| ) with spatial vari- 
ations that inc rease on larger angular sc ales, especially 
from 1° to 10° (|Adams fc SkrutskieifT99l) . We therefore 
conduct observations on a sounding rocket flight, at alti- 
tudes above the layers in the atmosphere responsible for 
airglow emission at characteristic altitudes of ~ 100 km. 



To measure the ~ 10' peak in the REBL spatial power 
spectrum, it is necessary to image an area of sky on the 
order of a square degree. While one can image a large 
field with a mosaic using a small field of view, this re- 
quires a highly stable instrument. A wide field of view 
allows a measurement using single exposures in the short 
time available on a sounding rocket flight. 

The REBL elec tromagnetic spectrum is predicted to 
peak at 1— 2(im (jCoorav et al.l I2004L iKashlinsk v et al.1 
12004 ) due to the redshift-integrated Lyman-a emission 
feature, with a decreasing spectrum at longer wave- 
lengths that depends on the history of reionization and 
the presence of free-free emission from ionized gas sur- 
rounding the first galaxies. Observations in the op- 
tical and near-IR should detect this spectrum, which 
is distinct from that of local foregrounds, namely ZL, 
stars, galaxies, scattered starlight (i.e. diffuse galactic 
light), and other Galactic emission. Though ideally the 
wavelength coverage would extend out to ~ 5 /zm, the 
key wavelengths for REBL science bracket the 1— 2 /im 
peak. Longer wavelength information can be obtained 
by cross-correlating CIBER data with overlapping wide- 
field Spitzer and AKARI maps. 

The local-galaxy fluctuations foreground is mitigated 
by masking galaxies down to a given flux threshold. The 
masking depth needed depends on the residual cluster- 
ing a nd Poisson fluctuatio ns of galaxies below the cutoff 
flux. lSullivan~et al. ( 2007j) measured galaxy clustering as 
a function of cutoff from a wide-field ground-based NIR 
survey catalog. We note that the REBL is best discrim- 
inated from low-redshift galaxy clustering and Poisson 
fluctuations at 10 arcminutes, as is evident in Figure 
[T] by comparing the REBL and galaxy clustering power 
spectra. Thus wide-field observations are also helpful for 
discriminating REBL from local galaxy fluctuations. 
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The flux cutoff needed to separate the optimistic REBL 
model from local galaxy fluctuations is ~ 17 th Vega mag- 
nitude at 1.6 /im, as is evident from the curves in Fig- 
ure Q] The spatial density of galaxies brighter than 
17 th Vega magnitude is A^(> S) = 500 galaxies per 
square degree. The cutoff required to remove galaxies 
well below the expected CIBER instrument sensitivity 
is ~ 23 rd Vega magnitude at 1.6 /xm, corresponding to 
N(> S) = 1.5 x 10 5 galaxies per square degree. Thus we 
find an angular resolution of 4—80 arcseconds is needed 
to remove galaxies in order to lose less than 25 % of the 
pixels from masking. 

Galaxy masking can be accomplished using ancillary 
observations with greater point source depth, masking 
pixels in the CIBER images below the CIBER point 
source sensitivity. The fields observed in the first two 
flights of CIBER, listed in Tablc[IJ allows source masking 
using deep companion catalogs obtained in ground based 
NIR observations. Detail s on first flight observa tions of 
these fields is available in iTsumura et al.l (|2010D . These 
fields have also been observed in a search for REBL fluc- 
tuations by AKARI and Spitzer at longer wavelengths, 
allowing for a cross-correlation analysis with CIBER. 

2. INSTRUMENT DESIGN 

The Imager instrument consists of two wide-field re- 
fracting NIR telescopes each with an 11cm aperture, 
combined with band-defining filters, a cold shutter, and a 
1024 x 1024 HgCdTe 2.5 /im Hawaii- 10 focal plane array. 
The Imager optics were designed and built by Genesia 
Corporation using the cryogeni c index of refraction mea- 
surements of lYamamuro et all ([20061 ) . A schematic of 
the assembly is shown in Figure [2] The assembly hous- 
ing the Imager optics are constructed from aluminum 
alloy 6061, and the lenses are made from anti-reflection 
coated Silica, S-FPL53 and S-TIL25 glass. The assembly 
is carefully designed to maintain optical alignment and 
focus through launch acceleration and vibration. The 
aluminum housing is hard black anodized to reduce re- 
flections inside the cryogenic insert and telescope assem- 
bly, with the exception of the static baffle at the front of 
the assembly which is gold plated on its external surface 
and Epner laser black coatecQ on its inner surface. This 
scheme serves to reduce the absorptivity of the baffle 
on the side facing warm components at the front of the 
payload section, and increase the absorptivity to NIR 
light on the inside. At the other end of the camera, a 
focal plane assembly is mounted to the back of the op- 
tical assembly and thermally isolated using Vespcl SP-1 
standoffs. The assembly includes a cold shutter and ac- 
tive thermal control for each detector. In addition, a 
calibration lamp system illuminates the focal plane in 
a repeatable way to provide a transfer standard during 
flight. The design of the calibration lamp system is com- 
mon to all of the CIB ER instruments and is presented in 
iZemcov et all (|2012D . 

The optical transmittance of the two Imager filters are 
shown in Figure [3] The filter stack is located behind the 
optical elements and in front of the focal plane assembly 
and cold shutter as shown in Figure [2] Each lens pro- 
vides additional filtering for wavelengths that are out of 

1 Manufactured by Teledyne Scientific & Imaging, LLC. 

2 This is a proprietary process of Epner Technology, Inc. 



band for both instruments, as their anti-reflection coat- 
ings transmit less than 1.5% of light with wavelengths 
shorter than 0.75 /im or longer than 2.0 /im. 

Table [2] summarizes the design properties of the op- 
tics and detector system, and the measured efficiencies, 
bands, and read noise for the two cameras. The optical 
efficiency is the product of the reflectance and absorption 
of the anti-reflection coated lenses taken from witness 
samples. The instrument performance is calculated in 
the appendix based on data from Table [2] and presented 
in Table |U 

Once assembled, the cameras mount to an optical 
bench shared with the LRS and NBS. The completed 
instrument section is then inserted into the experiment 
vacuum skin. Like the other CIBER instruments, the 
Imager optics are cooled to ~ 80 K to reduce their in- 
band emission using a liquid nitrogen cryostat system. 
IZemcov et al.l (j2012f i describes the various payload con- 
figurations used in calibration and in flight which allow 
both dark and optical testing in the laboratory. 

3. INSTRUMENT CHARACTERIZATION 

REBL fluctuation measurements place demanding re- 
quirements on the instrument, including the detector 
noise properties, linearity and transient response, opti- 
cal focus, control of stray radiation, and knowledge of 
the flat field response. We have carried out a series of 
laboratory measurements to characterize these proper- 
ties. 

3.1. Dark Current 

The detector dark current is measured in both flight 
and laboratory configurations by closing the cold shut- 
ters, which attenuate the optical signal by a measured 
factor of ~ 10 3 . Array data are acquired at 6.8 /xs per 
pixel sample, so that the full array is read in 1.78 s. The 
pixels are read non-destructively, and integrate charge 
until reset. The integration time may be selected, but the 
flight integrations are typically ~ 50 s. To maximize the 
signal-to-noise ratio, for each pixel we fit the measured 
output voltage to a slop e and an offset as described in 
iGarnett fc Forrest! (|1993l ). All CIBER Imager data arc 
analyzed using this method, except where noted. 

The measured dark current also depends on the de- 
tector thermal stability. For the Imagers we require 
dark current stability of 0.1e _ /s, which is equivalent to 
±100 /iK/s given a temperature coefficient of 1000 e~/K. 
The Imager detector arrays are controlled to ±10 /iK/s 
both in the lab and i n flight, exceeding this specification 
(|Zemcov et alj 120121 ). In the flight configuration with 
the cold shutter closed and the focal plane under active 
thermal control, we achieve ~ 0.3c _ /s mean dark cur- 
rent, as shown in Figured! The dark current is measured 
frequently before launch as a monitor of the instrument 
stability and is entirely consistent with the dark current 
measured in the laboratory. The stability of the dark 
current from run to run indicates the dominant contrib- 
utor to dark current is the array itself, as opposed to 
temperature or bias drift. 

3.2. Noise Performance 

Measuring the REBL spatial power spectrum requires 
a precise understanding of the noise properties of the 
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Fig. 2. — Schematic and photograph of the CIBER imaging camera. Light enters the optical system at left and is imaged to the focal 
plane at right. A fixed baffle is used to reduce scattering on the first optic. The I mager assembly emplo ys a fiber-fed calibration lamp 
system, band-defining and blocking filters, and a focal plane assembly as described in Zemcov et al. (2012). Both Imager assemblies used 
in CIBER are identical except for their band defining filters, set with AA/A ~ 0.5 bandpasses centered at 1.1/xm and 1.6 /im, roughly 
corresponding to astronomical I and H band. The photograph shows a fully assembled Imager in the lab. The entire assembly mounts to 
the CIBER optical bench when installed in the payload and operates at ~ 80 K. 



array. The array noise introduces a bias that must be 
accounted and removed in auto-correlation analysis, and 
determines the uncertainty in the measured power spec- 
trum. The instrument sensitivity shown in Figure Q] as- 
sumes the noise over the array is uncorrelated between 
pixels. Unfortunately, HgCdTe array s exhib it correlated 
noise, as described bv lMoselev etafl (I2010D . This noise 
is associated with pickup from the clock drivers to the 
signal lines, with 1// noise in the multiplexer readout, 
and depending on the implementation, with 1 // noise on 
the bias and reference voltages supplied to the array. 

3.2.1. Noise Model 

We characterized array noise using dark laboratory im- 
ages and data obtained just prior to flight. We first took 
a series of dark integrations to characterize the noise be- 
havior similar to the ^ 50 s integrations used in flight. In 
the left hand panels of Figure [5] we show the two dimen- 
sional power spectrum of the difference of two consec- 
utive 50 s laboratory integrations. The spectrum shows 



enhanced noise at low spatial frequencies along the read 
direction that is largely independent of the cross-read 
spatial frequency, symptomatic of correlated noise in the 
readout. 

We then generate an estimate of the noise by construct- 
ing time streams for the array readout. First, we deter- 
mine the best fit slope and offset for each pixel. We then 
subtract this estimate of the photo current signal in each 
pixel in each frame. Finally, we form a sequence of data 
for each of the four readout quadrants in the order that 
the readout addresses individual pixels. An example of 
time-ordered data and its noise spectrum is shown in 
Figure [6l ex hibiting excess noise b ehavior similar to that 
described in lMoselev et all (|2010l ). 

The correlated noise in the readout may reduce the in- 
flight sensitivity, and must be modeled to remove noise 
bias in the auto-correlation power spectra. While a full 
description of a noise model of the flight data is outside 
the scope of this paper, we can generate a model confined 
to the noise properties of the arrays observed in labora- 
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Fig. 3. — 1.1 fim and 1.6 fim Imager filter responses. These 
curves represent the transmission of the optical stack which in- 
cludes band defining and blocking filters as well as 5 anti-reflection 
coated lenses. This response does not include the response of the 
detector array, which typically cuts off at ~ 900 nm for a Hawaii- 1 
array with a sapphire substrate (Mark Farris, private communica- 
tion) . 



TABLE 2 
Imager Instrument Properties. 





1.1 fim Band 


1.6 fim Band 


Units 


Wavelength Range 


900-1320* 


1150-2040 


nm 


Pupil Diameter 


110 


110 


mm 


F# 


4.95 


4.95 




Focal Length 


545 


545 


mm 


Pixel Size 


7x7 


7x7 


arcsec 


Field of View 


2.0 x 2.0 


2.0 x 2.0 


deg 


Optics Efficiency 


0.90 


0.90 




Filter Efficiency 


0.92 


0.89 




Array QE 


0.51 


0.70 


** 


Total Efficiency 


0.42 


0.56 




Array Format 


1024 2 


1024 2 




Pixel Pitch 


18 


18 


fim 


Read Noise (CDS) 


10 


9 


c~ 


Frame Interval 


1.78 


1.78 


s 



* We assume a 900 nm cut-on wavelength from the 
Hawaii-1 substrate. 

** Array QE is estimated from QE measured at 2.2 fim 
for each array and scaled based on the response of 
a typical Hawaii-1. 



Fig. 4. — CIBER Imager dark currents for both cameras. The 
mean dark current is 0.3 e — /s, which is consistent with the man- 
ufacturer's specifications for Hawaii-1 arrays operating near LN2 
temperature. 

3.2.2. Estimated Flight Sensitivity 

To calculate the effect of correlated noise on the final 
science sensitivity, we take our sequence of dark labora- 
tory images, calculate the two dimensional power spec- 
trum, and apply a two-dimensional Fourier mask that re- 
moves modes sensitive to the excess low frequency noise. 
We remove these modes because they have a phase co- 
herence in real data that is not fully captured by the 
Gaussian noise model. After Fourier masking, we cal- 
culate the spatial power in logarithmic multipole bins. 
We then evaluate the standard deviation in the spatial 
power among eight dark images, and refer this to sky 
brightness units using the measured calibration factors 
in Table Because the laboratory data do not have 
appreciable photon noise, we add an estimate of uncor- 
rected photon noise from the flight photo currents. We 
compare this empirical determination of t he noise wit h 
the naive sensitivity calculation in Figure 171 (jKnoxll 1995| ) . 
The empirical noise is close to the naive calculation on 
small spatial frequencies, but is degraded by correlated 
noise on large spatial scales. However the instrument is 
still sufficiently sensitive to easily detect the optimistic 
REBL power spectrum. For future experiments, one may 
address the reference pixels in Hawaii-RG arrays to mit- 
igate the effects of correlated noise. 



tory testing. This model is generated by producing a 
Gaussian noise realization of the power spectrum given 
in Figure [6] This is used to generate random realizations 
of time ordered data. These data are mapped back into 
raw frames, and fit to slopes and offsets to determine 
the images for a full 50 s integration. To generate images 
like those shown in Figure [5j we generate multiple im- 
ages and display the difference of two 50 s images. This 
formalism will be extended to the flight data by adding 
photon shot noise from the astrophysical sky, and cor- 
recting for source masking, in a future publication. 



3.3. Detector Non-linearity and Saturation 

The Imager detectors have a dynamic range over which 
the response tracks the source brightness in a linear fash- 
ion. As is typical for Hawaii-1 detectors, the full well 
depth is measured to be ~ 10 5 e~; however, the detectors 
begin to deviate from linearity well before this. In order 
to flag detector non-linearity, we find pixels with differ- 
ent illumination levels and track their behavior during 
an integration. Figure [5] shows the typical response of a 
pixel to a bright ~ 3500 c~/s source over time. This plot 
shows a deviation from the linear model which is large 
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Fig. 5. — Images (top) and two dimensional power spectra (bottom) of the difference between two dark images, each obtained in a 50s 
integration. The upper and lower left hand panels show the image and power spectra of data take n m inutes before flight, while the two 
right hand panels show the same for random realizations using the noise model presented in Section I3.2I The spatial scale of these images 
has been restricted to 250 X 250 pixels to better show the spatial structure. In both cases the read direction is horizontal along pixel rows. 
The vertical structure in the two dimensional power spectra shows increased noise power in the read direction on scales > 50 pixels. The 
noise model accurately captures this behavior, both in real and Fourier space. 



at half the full well depth. Except for a few bright stars, 
Imager flight data are well within the linear regime. Pix- 
els with an integrated charge greater than 7000 e~ have a 
non-linearity ~ 1% are simply flagged and removed from 
further analysis, amounting to a pixel loss of < 0.5% over 
the array. 

3.4. Focus and Point Spread Function 

CIBER is focused in the laboratory by viewing an 
external collimated source through a vacuum window. 
Early on in focus testing we found that the best focus po- 
sition depended on the temperature of the optics. Ther- 
mal radiation incident on the cameras can heat the front 
of the optics and affect their optical performance due to 
both differential thermal expansion and the temperature- 
dependent refractive index of the lenses. We reduced 
the incident thermal radiation by installing two fused 
silica windows in front of the cameras for laboratory test- 
ing. The cold windows themselves are 125 mm diameter, 



5 mm thick Si02, operating at a temperature of 120 K, 
and have 1/10 surface flatness and < 5" wedge. As de- 
scribed in lZemcov et al.l (|2012[ ). these windows are ther- 
mally connected to the radiation shield to direct the ab- 
sorbed thermal power to the liquid nitrogen tank instead 
of routing the power through the optical bench where it 
would produce a temperature gradient across the optics. 

With the cold windows in place, we measure focus us- 
ing a collimator consisting of an off-axis reflecting tele- 
scope with a focal length of 1900 mm, a 235 mm unob- 
structed aperture, and an 8 /im pinhole placed at prime 
focus. Since the focus position of the instruments is fixed, 
we scan the pinhole through the focus position of the 
collimator to find the displacement from collimator best 
focus at which each Imager has its best focus. This pro- 
cedure is repeated at the center of the array, the corner 
of each quadrant, and in the center again as a check of 
consistency. Figure |H] shows data from such a test. If 
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Fig. 6. — The upper panel shows 30 ms of signal-subtracted time 
ordered data from the 1.6 /an Imager. The lower panel shows the 
noise spectrum derived from a longer such time series of reads 
over 50 s. The noise increases at ~ 10 kHz, visible in the time 
stream in the upper panel as the characteristic scale of the noise at 
~ 0.5 ms. The ringing visible in the power spectrum below 10 kHz 
corresponds to the harmonics of the clock signals used to address 
the array. 

the focal plane focal distance is found to be outside the 
±80 /im focal depth of the Imagers, we mechanically shim 
the focal plane assembly to the best focus position and 
remeasure the focus. We verify the focus position before 
and after prc-flight vibration testing, performed for each 
flight, to ensure that the focus will not change in flight. 

We measure the point spread function (PSF) in flight 
using stars as point sources. Given the large number of 
sources detected in each field, a measurement of the aver- 
age PSF across the array can be obtained by fitting all of 
the bright sources. In fact, because the astrometric solu- 
tion of the images allows us to determine source positions 
more accurately than a single pixel, and because the pix- 
els undersample the PSF of the optics, stacking sources 
gives a more accurate determination of the central PSF. 
To generate the stack, the region containing each source 
is re-gridded to be 3 x finer than the native resolution. 
The finer resolution image is not interpolated from the 
native image, rather, the nine pixels which correspond 
to a single native pixel all take on the same value. How- 
ever, when we stack the re-gridded point source images 
we center each image based on the known source posi- 
tions, and thus the stacked PSF is improved using this 
sub-pixel prior information. 

To measure the extended PSF, we combined data from 
bright sources, which saturate the PSF core, with faint 
sources that accurately measure the PSF core. We gen- 
erate the core PSF by stacking sources between 16.0 
and 16.1 Vega magn itudes from the 2MASS catalog 
(|Skrutskie et al.ll2006ft . which provides a set of sources 
that are safely in the linear regime of the detector. The 
source population is a combination of stars and galax- 
ies, however with 7" pixels, galaxies are unresolved. As 



a check, this same analysis was repeated for sources be- 
tween 15.0 and 15.1, and 17.0 and 17.1 Vega magnitudes. 
The PSF generated from these magnitude bands agreed 
with the nominal PSF. 

To measure the extended PSF, we stack bright sources 
between 7 and 9 Vega magnitudes from the 2MASS cat- 
alog. Since these bright sources are heavily saturated, 
the best fit Gaussian is only fit to the outer wings for 
normalization. After the core and extended PSFs are 
created, we find they agree well in the region between 
r <~ 13 arcsec, inside of which the bright sources are sat- 
urated, to r ~ 30 arcsec, where the faint sources are lim- 
ited by noise. 

We synthesize the full PSF by matching the amplitudes 
of the core and extended PSFs in the overlap region, 
producing the smooth two dimensional PSF shown in 
Figure[TUl The radial average of this full PSF is shown in 
Figure [TT] and highlights that the core PSF is consistent 
with the laboratory focus data. However, the extended 
PSF deviates significantly from this approximation and 
is better described by a Voigt profile shape, characteristic 
of scattering in the optical components. 

The extended PSF is essential for determining the ap- 
propriate mask to apply for bright sources. The diam- 
eter of the PSF mask is adjusted based on the bright- 
ness of the source, and pixels above a given flux are 
cut. The cut is calculated by simulating all sources in ei- 
ther the 2MASS or Spzteer-NDWFS catalogs using their 
known fluxes and the Imager PSF. The cut mask is gener- 
ated by finding all points on this simulation with fluxes 
> 3.3nW/m 2 /sr and 1.8nW/m 2 /sr at 1.1 and 1.6 fxm, 
respectively. This masking algorithm retains ~ 50 % of 
the pixels for a cutoff of 18 Vega mag, and ^30% of the 
pixels for a cutoff of 20 Vega mag. To test the cutoff 
threshold, we simulate an image of stars and galaxies 
and find that, cutting to 20 mag, the residual spatial 
power from masked sources is < 8 x 10~ 2 nW 2 m~ i sr~ 2 
at I = 10 4 , comparable to the instrument sensitivity 
shown in Figure [7] 

3.5. Off- axis response 

The Imagers must have negligible response to bright 
off-axis sources, including the ambient-temperature 
rocket ski n and shutter door, and the Earth. As de- 
scribed in lZemcov et al.l ([2012D . we added an extendable 
baffle to eliminate thermal emission from the rocket skin 
and experiment door, heated during ascent by air fric- 
tion, from illuminating the inside of the Imager baffle 
tube and scattering to the focal plane. 

We measured the off-axis respo nse of the fu l l baffl e sys- 
tem following the methodology in lBock et al.l (|1995D . We 
replaced the Hawaii- 1 focal plane array with a single op- 
tical photo diod^3 detector an d measured the respon se to 
a distant chopped source (see lTsumura et ahll2012l for a 
complete treatment of the measurement). The telescope 
gain function, 

g(9) = ^G{0), (1) 

where O is the solid angle of the detector and G(0) is 
the normalized response to a point source is the quan- 
tity of interest for immunity to off-axis sources in sur- 

3 Hamamatsu Si 10 X 10 mm 2 detector part number S10043. 
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Fig. 7. — The Imager sensitivity to REBL fluctuations. The left hand panel shows the estimated sensitivity for the 1.1 fim channel, and 
the right for the 1.6/^m channel. In addition to the curves taken from Figure[l] we show the sensitivity derived from laboratory data for 
both bands as described in the text using the same i binning as the naive sensitivity estimate shown by the orange curve. The black curve 
is an estimate of the flight sensitivity, combining measured laboratory noise from an ensemble of 50 s integrations, added with uncorrclated 
photon noise derived from the flight photo currents. This estimate is for a single 50 s integration, and does not include the effects of noise 
in the flat field or the loss of pixels from galaxy masking. 



120 
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Fig. 8. — Integrated signal as a function of time for a typical 
Imager pixel. The black data show subsequent reads of the Imager 
detector for an incident brightness of ~ 3500 c _ /s. The dashed red 
line shows the linear model matching the slope of the first 10 s of 
the integration. Finally, the blue line is a fit to the model from 
Bicsiadzinski ct al. (2011), which agrees well with the data. 

face brightness measurements (jPage et al.l 120031) and is 
independent of the optical field of view. The gain func- 
tion was measured for three baffle configurations and is 
shown in Figure 1121 The improvement from blackening 
the baffle tube and adding an extendable baffle section 
is notable for angles 9 > 20°. The stray light level from 
the Earth is given by 

W = ^ / 9{6)I®{6,4>)dSl, (2) 

where I§ is the surface brightness of the Earth, and 
Astray is the apparent surface brightness of stray light re- 
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Fig. 9. — The variation of the PSF width measured in the labo- 
ratory as a function of collimator focus position Ax shifted away 
from its best focus position. At each collimator position we mea- 
sure the PSF by fitting a Gaussian and determining its full width 
at half maximum (FWHM) and uncertainty. The points show the 
data and the black line the best fit parabola to the points, yielding 
the best estimate of the focus position of the Imager instrument. 
The curve is consistent with the //4.95 focal ratio, where the array 
pixels are 18 X 18 fim and subtend 7x7 arcseconds on the sky. 

ferred to the sky. Fol lowing the calculation described in 
iTsumura et al.1 ()2012l ). we estimate that during the sec- 
ond flight CIBER observations of the fields listed in Table 
[TJ where the Earth's limb is > 72° off-axis, the stray light 
level is calculated to be 2 nW/m 2 /sr and 1 nW/m /sr in 
the 1.6 and 1.1 /im channels, respectively. 

This level of stray light is quite small but not com- 
pletely negligible, and potentially problematic in an 
anisotropy measurement depending on its morphology 
over the field of view. To quantify how stray light af- 
fects our measurements, we calculated the spatial power 
spectrum of the difference between two images, BootesA 
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Fig. 10.— The 1.1 /jm (left) and 1.6 /im (right) Imager PSFs 
measured using stacked flight images from a combination of bright 
and faint sources as described in the text. The Imager PSF has 
a bright core with a faint extension to r ~ 1', and is circularly 
symmetric. 

- BootcsB which are separated by only 2° on the sky and 
taken at nearly the same Earth limb avoidance angle, and 
Bootes A - NEP, from second flight data (see section [5]) . 
We find that the power spectra of these differences are 
the same to within statistical noise, and that the spatial 
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Fig. 11. — The radial profile of the 1.1 fj,m (left) and 1.6 /xm 
(right) Imager flight PSFs from Figure llOl (blue circles). The red 
curve shows the best fit Gaussian to the PSF core, while the black 
curve shows a best fit Voigt (i.e. the convolution of a Gaussian 
and Lorentzian) function to the extended PSF. This is indicative 
of scattering in the optical components. Finally, the black dash- 
dotted line shows the HWHM of the PSF, which matches the value 
measured in the laboratory. 

fluctuations of the stray light signal are negligible. 

We plan to observe these fields again in future flights 
at different Earth limb avoidance angles, including angles 
greater than 90°. The cross-correlation of such images 
from different flights is highly immune to residual stray 
light. 

3.6. Flat Field Response 

The instrumental flat field, which is the relative re- 
sponse of each detector pixel to a uniform illumination 
at the telescope aperture, is determined in flight by av- 
eraging observations of independent fields. Additionally, 
the flat field can be independently measured in the lab- 
oratory before and after flight as a check for systematic 
error. The laboratory flat field response is measured by 
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Fig. 12. — The Imager telescope gain function, measured with 
the anodized fixed black baffle tube used in the first flight (dotted 
black line), an improved fixed baffle tube with a better laser black 
optical coating (Epncr Technology Inc., dashed blue line), and the 
combination of the improved fixed baffle with an extendable baffle 
used in the second flig ht (solid red line). D etails of the optical 
baffling can be found inlZcmcov ct al. (20j3). 

illuminating the full aperture of a camera with the out- 
put of an integrating sphere. The sphere is illuminated 
with a quartz-tungsten halogen lamp which is filtered to 
produce an approximately solar spectrum at the output 
of the sphere, mimicking the spectrum of ZL. 

The sphere was measured by the manufacturer to have 
uniformity as a function of angle to better than 5 x 1CP 3 
over 10° x 10°. We scanned a small collimathig telescope 
with a single pixel over the aperture, and determined that 
the sphere has angular uniformity to better than 1 x 10 -3 
over the 2° x 2° Imager field of view. We also measured 
the spatial uniformity over the output port and saw no 
evidence of non-uniformity to < 7 x 10 -3 over an 11 cm 
aperture. 

To eliminate any effects from vacuum and thermal win- 
dows, we house the integrating sphere inside a vacuum 
chamber which mates t o the front of the cryo stat in place 
of the shutter door (see lZemcov et al.l (|2012[ ) for details). 
Light is fed into the sphere from outside of the vacuum 
box so that the lamp can be chopped at the source, allow- 
ing us to remove the thermal background. An example 
flat field measurement for the 1 . 1 /xm camera is shown in 
Figure [T3J 

The laboratory data arc fitted over a limited period 
of the integration following array reset so as to avoid an 
appreciable error from non-linearity, as described in sec- 
tion [331 taking into account the minimum well depth of 
all pixels in the array. The instruments have a residual 
response to thermal infrared radiation in the laboratory 
with a typical photo current of 600e~/s in the 1.6 //m 
array, which therefore limits the linear integration pe- 
riod to ~ 5 s. We obtained interleaved data with the 
source on and off to monitor and subtract this thermal 
background. After accounting for these effects, the final 
statistical accuracy of the laboratory flat field images 
shown in Figure [TBI is 1.6% per pixel. Laboratory flat 
fields were measured before and after the second flight 
to quantify the reproducibility of the lab flat field re- 
sponse. We binned 1.6 /xm camera laboratory flat field 



images into 64 (15 x 15) arcminute square patches in 
order to reduce statistical noise, and found the binned 
images agree to < l%(lcr). The agreement between the 
flight and laboratory flat fields requires a full reduction 
of the flight data and will be presented in a future science 
paper. 

4. MODIFICATIONS FOLLOWING THE FIRST FLIGHT 

The Imagers were flown on the CIBER instrument on 
a Terrier Black Brant sounding rocket flight from White 
Sands Missile Range in 2009 February. Many aspects of 
the experiment worked well, including the focus, arrays 
and readout electronics, shutters, and calibration lamps. 
However, we also found several anomalies that led to 
modifications for subsequent flights. 

4.1. Thermal Emission from the Rocket Skin 

The instruments showed an elevated photon level dur- 
ing the flight due to thermal emission from the rocket 
skin, heated by air friction during ascent, scattering into 
the optics. The edge of the skin near the shutter door 
can directly view the first optic and the inside of the 
static baffle. This thermal response was pronounced at 
long wavelengths, as traced by the LRS (jTsumura et al.l 
l2010l l. The 1.6 (im Imager was more affected by thermal 
emission than the 1.1 /xm Imager, as expected from its 
longer wavelength response, giving 40 and 7 times the 
predicted photo current, respectively. 

The measured thermal spectrum with the LRS should 
not produce a significant photo-current in the 1.1 /xm 
Imager, as the band is supposed to cut off at 1.32 /im. 
The excess photo-current indicates the 1.1 /xm Imager has 
some long wavelength response. The array response may 
continue somewhat beyond 2.5 /xm, as the band-defining 
filters provided blocking out to just 2.5 /xm and th en open 
up. Also as with the NBS (jKorngut et al.l 120121 ) the fil- 
ters may not attenuate scattered light at large incident 
angles as effectively as at normal incidence. The bright- 
ness observed by the 1.6 /xm Imager is 6 times higher 
than the band-averaged LRS brightness. This could be 
due to a combination of the higher stray light response 
in the 1.6 /xm Imager, and the filter blocking issues men- 
tioned above. We installed an additional blocking filter 
providing < 0.1% transmittance from 2.4 /im to 3.0 /xm 
for both imagers. 

We modified the front of the experiment section to bet- 
ter control the thermal and radiative environment at the 
telescope apertures. Most notably, we added extendable 
baffles to each of the instruments to eliminate all lines 
of sight from the sk in to the optics or the inside surfaces 
of the baffle tubes. iZemcov et al.l (|2012t) details the de- 
sign of these baffles and the other changes made to the 
experiment section front end. Thermal emission is not 
detectable in the Imagers in the second flight, and is at 
least 100 times smaller than the first flight in the LRS 
data. 

4.2. Rings and Ghosts from Bright Sources 

During analysis of the first flight data, we discovered 
that bright objects outside of the Imager field of view cre- 
ate diffuse rings in the final images, as shown in Figure 
[T4l Upon further analysis, we found that each of these 
rings was centered on a bright star outside the geometric 
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Fig. 13. — The 1.1 /an and 1.6 /xm Imager flat fields as measured in the lab using the apparatus described in Zcmcov ct al. (2012j). The 
average response has been scaled to 1.0 in this image, which shows the typical relative responsivity performance of the Hawaii-1 arrays in 
conjunction with the optics. The RMS variation in the pixel responsivities is 0.09 at 1.1 /an and 0.12 at 1.6 fim. 



field of view. The rings were caused by reflections off in- 
ternal elements of the telescope assembly, as illustrated 
in Figure 1151 There are two general classes of rings in 
the first flight images, though the second class contains 
two distinct populations; we denote these ring popula- 
tions 1, 2 and 3 below. Table [3] gives details of the ring 
populations including their angular extent and coupling 
coefficients. 

Population 1 rings are generated by reflections off a 
lens mounting flange (Figure 1 1 5|) , and are produced by 
bright sources between 3.4° and 6.6° off-axis. These rings 
also have the strongest optical coupling, with an inte- 
grated flux in the ring a few tenths of percent of the 
incident source flux. Given their large acceptance angle, 
stars brighter than 4 th magnitude are sufficiently abun- 
dant to generate multiple bright rings. 

Following their discovery in the first flight data, we 
measured the population 1 rings and searched for other 
optical reflections in the laboratory. We illuminated each 
Imager aperture with collimated light and then scanned 
the angle of incidence of the collimated beam up to 25° 
off-axis. The first set of measurements confirmed the 
existence of the population 1 rings, and allowed the dis- 
covery of the second class of fainter rings. 

The second class of rings is comprised of two sub- 
populations which are both generated by reflections off 
the lens tube and lens support fixtures at the front of 
the optics assembly (Figure IT5|) . These rings have flux 
coupling coefficients similar to, but slightly less than, the 
population 1 rings, but have much larger solid angles on 
the array and so produce smaller per pixel brightness. 
Together, population 2 and 3 rings are caused by bright 
sources 6.7° to 13.2° off-axis. These rings are not readily 
visible in the images from the first flight, though their 
presence was verified in the lab after flight. 
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Fig. 14. — 1.1 (j,m image of the Bootes A field from CIBER's 
first flight showing rings which were later traced to reflections off 
components inside the Imagers, namely the lens mounts and in- 
strument walls. As a guide the brightest rings are indicated with 
arrows. There are three separate populations of reflections which 
produce these rings. All sources which fall into their angular re- 
sponse regions will produce a ring, though only sources brighter 
than magnitude ~ 4 produce rings which are visible by eye. These 
rings produce excess power in the science power spectrum, but were 
eliminated by modifying the optics for the second flight. 



Given the acceptance angles, star number counts and 
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TABLE 3 

First Flight Imager Ring Parameters. 

Ring Type / during (<!>)/ J I 

Pre-fix Post-fix (3c) Reduction in Ci 

1.1/Ltm Imager 

1 3A° 6S° 2.2 x 1CT 3 < 2.6 x 10~ 6 > 7 x 10 5 

2 6.7° 8.8° 2.7 x 1(T 4 < 1.5 X 10 -6 > 3 x 10 4 

3 11.2° 13.2° 6.6 x 10~ 4 < 1.6 x 10~ 6 > 1 x 10 5 

1.6/xm Imager 

1 3A° 6S° 4.1 x 10" 3 < 3.0 x 10~ 6 > 1 x 10 s 

2 6.7° 8.8° 3.5 x 1CT 4 < 1.9 x 1CT 6 > 3 x 10 4 

3 11.2° 13.2° 1.3 x 1CT 3 < 4.3 x 1CT 6 > 9 x 10 4 




Fig. 15. — Ray trace from an off-axis source which produces the 
rings observed at the focal plane. The first class of rings (labeled 
as 5° in the Figure) are caused by glancing reflections off a flange 
supporting the back lens. The second class of rings (labeled as 7.5°) 
is produced by glancing reflections off flanges and lens holders in 
the front set of optics. For the second flight, these surfaces were cut 
back and grooved to reduce the glancing reflectance, removing the 
rings to a negligible level, as verified by laboratory measurements. 

the quality of the ancillary data, the first set of rings are 
sufficiently bright to be modeled and masked from the 
first flight images. However, the second set of rings have a 
more complex morphology and fainter surface brightness, 
and are more difficult for us to confidently account for in 
the images. 

To understand the systematic error associated with the 
population 2 and 3 rings, we modeled their effect by con- 
volving the measured laboratory response with an off- 
axis star catalog for each field, and calculated the spatial 
power spectrum of the resulting images. These rings, 
if left unmasked, produce power above the instrument 
sensitivity level, as shown in Figure [T6l 

To remove the rings entirely, we made the optical sim- 
ulation shown in Figure 1151 Following characterization 
of the rings, the Imager optical assemblies were disas- 
sembled. The components responsible for the rings were 
grooved or cut back and re-anodized. The Imager optics 
were then reassembled, and the off-axis measurements 
were repeated. We did not observe any rings following 



these modifications. We place upper limits on the ring 
coupling factors shown in Table [3] which are based on 
the uncertainty in the integrated surface brightness over 
the nominal ring solid angles from the laboratory mea- 
surements. We propagated these upper limits through 
the model to produce synthetic images and then power 
spectra. The estimated reduction in the power spectrum 
from the first class of rings are given in Table [3l We 
find that the effect on the power spectrum is negligible 
compared with the instrument sensitivity after the optics 
modifications. 

5. INSTRUMENT PERFORMANCE FROM THE SECOND 
FLIGHT 

The Imagers were flown on the CIBER instrument on 
a second sounding rocket flight in 2010 July. All aspects 
of the experiment performed well. We found no evidence 
of bright thermal emission from the rocket skin in cither 
of the Imagers. We did not observe rings in the flight 
images. While the science data are still being analyzed, 
we summarize the observed brightness and array photo- 
currents in Table |U Unfortunately, it is difficult to esti- 
mate the full in-flight sensitivity in the power spectrum 
without a noise estimator that accounts for correlated 
noise in the presence of sources and masking. Therefore 
we estimate the in-flight per-pixel sensitivities by evalu- 
ating the noise in the flight difference images (see Section 
13. 2[) . The corresponding per pixel surface brightness sen- 
sitivities, and point source sensitivities using a 2 x 2 pixel 
aperture, are listed in Table |4j Our estimated sensitiv- 
ity to the spatial power spectrum is shown in Figure [7] 
based on the variance of the power spectra of an ensemble 
of dark laboratory images combined with flight photon 
noise. 

We scale the photo currents in Table [4] to sky bright- 
ness units using a calibration based on point sources ob- 
served in flight. We stacked sources with flux between 
16.0 and 16.1 Vega magnitudes in the 2MASS catalog and 
integrated over the stacked image to account for the ex- 
tended PSF. We converted this point source calibration 
to surface brightness using the pixel solid angle, giving 
the calibration factors in Table [4] 

6. CONCLUSIONS 

We have designed and tested an imaging instrument 
optimized to search for the predicted spatial and spectral 
signatures of fluctuations from the epoch of reionization. 
The instrument demonstrates the sensitivity needed to 
detect, or place interesting limits upon, REBL fluctua- 
tions in the short observing time available in a sounding 
rocket flight. We have carried out a comprehensive lab- 
oratory characterization program to confirm the focus, 
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Fig. 16. — Simulated power spectra for the second class of rings for both Imager instruments, 1.1 (im (left) and 1.6 (im (right). These 
spectra were computed given the ring parameters in Tablcl3land the known star fluxes and positions near the CIBER fields. The instrument 
sensitivity is the the same as modeled in Figure [T] The amplitude of the power spectrum of the rings is different for each field because 
of the differing stellar populations near each, but similar between the bands because of the typical color of stars. For the second flight, 
the level of ring contamination is well below the instrument sensitivity, based on upper limits obtained in the laboratory following the 
modifications to the optics described in the text. The upper limit is shown for SWIRE, the most demanding field. 



TABLE 4 

Calculated and Second Flight Sensitivities in a 50 s Observation. 





1 . 1 fj,m 


Imager 


1.6 fim Imager 






Predicted 


Achieved 


Predicted 


Achieved 




Sky brightness 


450 


420 


300 


370 


nW m~ 2 sr _1 


Photo current 


4.4 


4.9 


8.2 


11.0 


e _ /s 


Responsivity 


10 


11 


28 


31 


me - s — 1 / nW m~ 2 sr _1 


Current Noise 


0.31 


0.35 


0.41 


0.45 


e - s — 1 (lcr/pix) 


6XI X 


31.7 


33.1 


15.1 


17.5 


nW m -2 sr _1 (lcr/pix) 


SF V 


18.5 


18.4 


18.2 


17.8 


Vega Mag (3cr) 



characterize the flat field response, perform an end-to- 
end calibration, and measure the stray light response and 
detailed noise properties. After a first sounding rocket 
flight in 2009 February, we modified the instrument to 
eliminate response to thermal radiation from ambient 
portions of the payload, and to reduce stray light to 
bright stars outside of the field of view. Scientific data 
from the second flight in 2010 July are currently under 
analysis, and the instrument demonstrated sensitivity 
close to design expectations. The instrument character- 
ization shows that systematic errors from the extended 
PSF, stray light, and correlated noise over the array are 
controlled sufficiently to allow a deep search for REBL 
spatial fluctuations. We recently completed a third flight 
in 2012 March that allows us to cross-correlate images at 
different seasons to directly assess any ZL fluctuations. 
The flight and recovery were successful, and a fourth 
flight is now planned. A successor instrument, with 3 or 
more simultaneous spectral bands and with higher sensi- 
tivity using a 30 cm telescope and improved Hawaii-2RG 
arrays, is currently in development. 

7. APPENDIX 

The calculated sensitivities in Table 21 Figure Q] and 
Figure [7] are based on a 50 s integration with the instru- 
ment parameters given in Table [2] The estimated photo 



current i p hot given by: 
iphot — XI\ 



/ r]Afl AX 
\ hv X 



[e"/s], 



(3) 



where AVL is the pixel throughput, 7/ is the total efficiency, 
XI\ is the sky intensity, and AA is the integral band- 
width. The term in brackets in Equation [3] gives the 
surface brightness calibration from e^/s to nW/m 2 /sr. 
The current noise over an integration with continuous 
sampling is given by: 



^phot 



-phot 



6Tn 



-•3cdsYT [e /s], (4) 



where <5Qcds is the correlated double sample read noise, 
T — 50 s is the integration time, and the frame rate To = 
1.78 s. The surface brightness sensitivity is therefore: 

6Xh = Siphat M^kxji [nW m ~ 2 sr ~ 1] - (5) 

Finally, the point source sensitivity is given by: 



5XF X = ^ phot ^S^ [nWm- 2 ], (6) 

where N P i X is the effective number of pixels that must 
be combined to detect a point source, and we have as- 
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sumed N p i X = 4. These per-pixel sensitivities are used 
to estimate the sensitivity on the power sp ectrum in Fig- 
ure [1] and Figure [7] using the formalism in iCoorav et al.l 
(|2004| ). The calculation assumes the noise in each pixel 
is independent, and ignores errors from source removal 
and flat-field estimation. 
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